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^ : ABSTRACT 

^ I A systematic study of the linear thermal stability of a medium subject to 

l/^ \ cooling, self-gravity and thermal conduction is carried out for the case when the 

^ I unperturbed state is subject to global cooling and expansion. A general, recursive 

O ! WKB solution for the perturbation problem is obtained which can be applied to 

^ I a large variety of situations in which there is a separation of time-scales for the 

"Th I different physical processes. Solutions are explicitly given and discussed for the 

case when sound propagation and/or self-gravity are the fastest processes, with 
O I cooling, expansion and thermal conduction operating on slower time-scales. A 

^ [ brief discussion is also added for the solutions in the cases in which cooling or 

I conduction operate on the fastest time-scale. The general WKB solution obtained 

. ^ I in this paper permits solving the problem of the effect of thermal conduction and 

^ ' self-gravity on the thermal stability of a globally cooling and expanding medium. 

I As a result of the analysis, the critical wavelength (often called Field length) above 

which cooling makes the perturbations unstable against the action of thermal 
conduction is generalized to the case of an unperturbed background with net 
cooling. As an astrophysical application, the generalized Field length is calculated 
for a hot (10^—10® K), optically thin medium (as pertains, for instance, for the hot 
interstellar medium of SNRs or superbubbles) using a realistic cooling function 
and including a weak magnetic field. The stability domains are compared with 
the predictions made on the basis of models for which the background is in 
thermal equilibrium. The instability domain of the sound waves, in particular, 
is seen to be much larger in the case with net global cooling. 

Subject headings: instabilities — intergalactic medium — ISM: clouds — meth- 
ods: analytical — waves — X-rays: diffuse background 
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1. INTRODUCTION 

Thermal instability often appears in astrophysical media which are undergoing a global 
process of cooling, expansion or contraction. A non- exhaustive list of examples includes: 
thermal instability in a cooling flow in a galaxy cluster (Mathews & Bregman 1978; Balbus 
& Soker 1989; Balbus 1991), in a stellar wind (Balbus & Soker 1989), or in the hot ISM 
(Begelman & McKee 1990), and structure formation in a collapsing protogalaxy (Fall & 
Rees 1985). In all those cases, the unperturbed state is not one of thermal or dynamical 
equilibrium, so that the stability and the time evolution of the perturbations cannot be 
studied by means of a classical normal-mode Fourier analysis in time. The situation is 
complicated by other physical processes which can be important for the unperturbed state 
or for the perturbation itself, like self-gravity, thermal conduction, and magnetic induction 
and diffusion. The stability criteria and the evolution of the perturbations directly depend on 
the relative time-scale of the latter processes with respect to each other and to the time-scales 
of cooling and expansion of the background in which they are developing. 

The astrophysical literature abounds in examples of thermal stability analyses and their 
application to specific problems, as in Parker (1953), Weymann (1960), Field (1965), Gold- 
smith (1970), Defouw (1970), Heyvaerts (1974), Mathews & Bregman (1978), Flannery 
& Press (1979), Balbus (1986, 1991), Malagoh, Rosner & Bodo (1987), Balbus & Soker 
(1989), Begelman & McKee (1990), Loewenstein (1990), Chun & Rosner (1993), Burkert 
& Lin (2000). A milestone was the paper by Field (1965). He considered a non-expanding 
homogeneous medium in thermal equilibrium (i.e., with cooling exactly compensating heat- 
ing), including thermal conduction and magnetic fields. The small-perturbation problem 
is then amenable to a normal-mode Fourier analysis in space and time. Additionally to 
magnetosonic waves, which can be strongly modified by the thermal processes, he found a 
non-oscillatory mode (called condensation mode) which may have an isobaric or isochoric 
character depending on the amount of cooling or conduction possible during one crossing 
time of a magnetosonic wave across the perturbation. Conduction becomes important at 
small spatial scales and, in general, has a strongly stabilizing effect. Another important 
step forward was taken by Mathews & Bregman (1978) and Balbus (1986), who consid- 
ered a time-dependent unperturbed medium. By using the entropy equation, they obtained 
stability criteria for the condensation mode in a medium with global cooling, not taking 
into account conduction, self- gravity or any background stratification. Balbus (1986), in 
particular, carried out a WKB analysis in space for radial perturbations in an expanding or 
contracting medium with spherical symmetry. Balbus & Soker (1989) relaxed the condition 
of homogeneous background and obtained solutions for a stratified medium. Balbus (1991) 
also included a background magnetic field. 
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The topic of thermal instability in astrophysical media seems to have been gathering 
momentum in the past several years. In spite of the foregoing, there is to date no systematic 
study of the stability of a magnetized medium subject to cooling, self-gravity and thermal 
conduction when the unperturbed state itself is undergoing expansion and cooling. The 
purpose of the present paper is to carry out such a study using, wherever possible, a WKB 
expansion in time. The treatment introduced here permits recovering many results of the 
previous literature in a unified manner and allows obtaining new ones. The non-magnetic 
and non-stratified case is considered in the present paper, which deals with the problem 
of a uniform (but time-dependent) background. However, the results of this paper can 
also be apphed locally to a weakly non- homogeneous background (see sect. 2.5). Several 
characteristic time-scales have to be considered, to wit, those associated with the background 
expansion, self-gravity, cooling, sound-wave propagation and conduction. When background 
cooling and expansion are not the fastest processes in the system, a WKB analysis in time 
can be carried out. General WKB solutions up to arbitrary order are then obtained and 
applied to the particular cases when the fastest time-scales are those associated with sound 
wave propagation, with sound and self-gravity simultaneously, or with conduction. 

New results obtained in the present paper are as follows. First, we obtain an expression 
for the critical length separating stable from unstable domains which generalizes the classical 
Field length to a system which, in the unperturbed state, is cooling (or otherwise) and 
expanding or contracting globally. Field (1965) showed the opposite effects of conduction 
and cooling on the stability of a system in static and thermal equilibrium in the case when 
the cooling function, £, fulfills [djC/dT]p < 0, and found that below a critical wavelength 
(the classical Field length) conduction stabilizes the system. Our study permits relaxing 
the condition of thermal equilibrium in the unperturbed state and including in the latter a 
global expansion or contraction. We may note that, for lack of a more accurate expression, 
different authors in the past have tentatively used the value calculated by Field even though 
their problem involved a medium undergoing global cooling (e.g., David & Bregman 1989; 
Pistinner & Shaviv 1995). In section 5, we calculate the generahzed Field length for a 
medium undergoing a realistic net cooling, and compare it with the classical Field length of a 
medium in thermal equilibrium with the same cooling but with a balancing constant heating. 
The results of the present paper also generalize those of Balbus (1986), who considered a 
medium undergoing global cooling in the unperturbed state, but ignored thermal conduction. 

Second, there is in the literature no simultaneous study of the thermal and gravitational 
linear stability of a medium undergoing net cooling. The combination of both phenomena 
is important, for example, for the structure formation in a protogalaxy (Fall & Rees 1985), 
where condensations previously formed by thermal instability in a cooling medium become 
gravitationally unstable. In section 3.2 we calculate how the acoustic-gravity modes of Jeans 
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(1902) are modified through conduction and coohng when the background itself is coohng 
globally. We also show the existence of a hitherto unknown condensation mode in this 
case. In all these solutions the perturbations can go through different stability regimes: in a 
cooling medium the generalized Field length and the Jeans length change with time whereas 
the wavelength of the perturbation remains constant. 

There is an alternative to a WKB expansion when trying to solve the problem tackled 
in this paper. Burkert & Lin (2000) have used a Taylor expansion in time for the reduced 
problem when there is no thermal conduction, self- gravity or background expansion. Yet, 
a Taylor expansion has the large disadvantage of being valid only at times small compared 
with the characteristic times of evolution of the background. This is in contrast to the 
WKB approach taken in the present paper, which is valid for arbitrary times as long as the 
perturbation does not grow beyond the linear regime. 

The present paper is organized as follows: in section 2, the equations of the problem 
are presented and the relevant time-scales introduced both for the background and the 
perturbations. Section 2.5 contains a discussion of the application of the WKB method 
to the present problem and provides a general WKB solution to the linear perturbation 
equations. In section 3, a collection of solutions for intermediate wavelengths is presented. 
The word intermediate is used in the sense that the wavelength is not small nor large enough 
for conduction or cooling, respectively, to be the dominant processes in the evolution. The 
case in which sound propagation is the fastest process is considered in section 3.1 whereas 
in section 3.2 both self-gravity and sound arc the fastest processes. Section 4 is devoted 
to the large and small wavelength regimes and, in particular, to the cases when cooling 
or conduction operate on the fastest time-scale. In section 5, the generalized Field length 
and the instability domain are computed for an optically thin medium at high temperature 
threaded by a magnetic field without dynamical effects (e.g., the hot ISM). The final section 
(sect. 6) contains a summary and further discussion of the results of the paper. 

2. THE EQUATIONS OF THE PROBLEM AND THEIR WKB SOLUTION 

This section discusses the physical processes relevant for the present problem, the gov- 
erning system of equations, and the method used to solve this system. The equations are 
given first in general (sec. 2.1) and then specialized for the homogeneous unperturbed medium 
(sec. 2.2) and for small perturbations to that medium (sec. 2.3). There is a considerable num- 
ber of different physical processes of importance for the present problem: the corresponding 
time-scales are listed and explained in section 2.4. A general WKB solution to the linear 
perturbation system is obtained and discussed in section 2.5. 
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2.1. Equations 

We consider an optically thin polytropic ideal gas. The fluid equations that describe 
the evolution of the system including thermal conduction and self-gravity are: 

| + pV-v = 0, (1) 

P^ + Vp + pV$ = 0, (2) 

_^,^+,V.v + p£-V.(xVT) = 0, (3) 

where d/dt = d/dt + \ - V is the Lagrangian time derivative, 7 is the polytropic index 
of the ideal gas, /i is the mean atomic mass per particle, TZ is the universal gas constant, 
C{p,T) (erg s~^ gr~^) is the net cooling function (including heating), x(p, T) is the thermal 
conduction coefficient and $ is the gravitational potential caused by the mass distribution 
of the system A$ = AnGp, with G the universal constant of gravitation. All other symbols 
have their usual meaning in hydrodynamics. Equations (1), (2) and (3) are the continuity, 
momentum and energy equations, respectively. The temperature is related to the other 
quantities by the ideal gas equation: p = {TZ/fi)pT. The change of /i due to recombination 
or ionization processes is not considered in this paper. This is fully acceptable for high 
enough temperature (T > lO^K), in which the medium is almost completely ionized, or in 
situations at lower temperature when the change in p is less important than the contribution 
of cooling and conduction. The thermal instability of a medium in thermal equilibrium with 
variable p has been considered by Goldsmith (1970) and Defouw (1970). Finally, note that 
the radiative flux in the diffusion limit (optically thick medium) is identical to the classical 
thermal conduction flux except for the associated conductivity coefficient (e.g., Mihalas & 
Mihalas 1984). Therefore, our treatment also allows for the use of radiative transport in 
the diffusion limit. 



2.2. Background Evolution 

As a basis for the small-perturbation analysis, we assume a homogeneous background 
which is expanding uniformly and cooling (or heating) by radiation, so that the unperturbed 
quantities only depend on time. The background quantities will be denoted with the sub- 
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script 0. The expansion is given by 

r = a(i)x, (4) 

where r is the Eulerian coordinate, x is the Lagrangian coordinate and a{t) is the expansion 
parameter. Using equation (4), the unperturbed velocity field is given by 



Vo(r,t) = -r, (5) 

where the dot over a symbol (as in d) indicates its time derivative. 
Equations (1) and (3) reduce for the background to 

-rr = -3-po , 6) 
at a 

|lo.T„^-(,-I)(3^.^|), P) 

respectively, where Cq — >C[po(^), T'o(i)]. The time dependence of a{t) must be specified from 
equation (2). The latter yields 

a — const (8) 

when self-gravity is negligible and 

4 

a = --TrGpoO (9) 

otherwise. Additionally, more general forms for a{t) are allowed for the local application of 
the perturbative analysis around a Lagrangian fiuid element in a non-uniform background, 
as will be explained at the end of section 2.5. 



2.3. Lineeirized Equations for the Perturbations 

To obtain a linearized system of equations for the perturbations we proceed as follows. 
First, we split each variable into unperturbed and perturbed components, indicating the 
latter with a subscript 1, like, e.g. 

Pi(r,t) p{r,t)-p,{t), (10) 
v,(r,t) v(r,t)-Vo(r,i), (11) 

and similarly for pressure and temperature. Eulerian linear perturbations are then taken in 
equations (1) through (3). Thereafter, the Eulerian divergence operator is applied to the 
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momentum equation and all equations are rewritten in terms of the Lagrangian coordinate 
X defined in equation (4). The resulting linear system has coefficients which depend on t but 
not on X. We then carry out a spatial Fourier analysis, with Fourier components proportional 
to exp(ik • x), so that k is the Lagrangian wave vector: 



pik(x, t) = pik(t) exp(i k • x) , 
vik(x, t) = vik(t) exp(?: k • x) , 



(12) 
(13) 



and expressions similar to (12) for all other scalar variables (pressure, temperature, etc). 
For simplicity, we eliminate the index k from the Fourier amplitudes in the following. The 
perturbed velocity is then split into its components parallel and normal to k: 



(14) 



Finally, the resulting linear system is simplified by repeated use of the background equations 
and elementary thermodynamic relations. We obtain: 



dt \po 



— = -ikvi , 



dvi 
'dt 



— Vi - I h — , 

a 1 Pq k Po 



d^fpi 
dt \pq 



^ ^ Po ^ po 



(15) 

(16) 
(17) 



with all coefficients of the perturbations being evaluated for the background values. As a 
result of the linearization, the time derivative d/dt in equations (15) through (17) is now given 
by d/dt = d/dt + vq-V, i.e., it is the material derivative as one moves with the unperturbed 
fluid. The other symbols have the following value: 



CP ^ (7-1)^^0 





Us = 
















(f)] 


p 



kcg ; 



AnGpo ■ 



(^_1)^^P 



7^po 



and 



CP - (7-1)^^0 







dT 


m 



where 



Po 
7— 
Po 



1/2 



and 



k 
a 



(18) 
(19) 

(20) 
(21) 
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are, respectively, the adiabatic background sound speed and the Eulerian wavenumber. The 
coefficients (18)-(20) provide characteristic times for our problem: their meaning and rele- 
vance will be discussed in the next section. One also obtains an independent equation for 
the perpendicular component of the velocity, viz.: 



i.e., the integral of the perturbed velocity along a closed material curve remains constant in 
time. 

The system (15)-(17) is a hnear ODE system in t for the vector of unknowns (pi/po, vi, pi/po)- 
Considering the ratio perturbation/background is the natural choice because the unper- 
turbed quantities depend on time. We will say that there is instability if either \pi/ Po\ or 
\Ti/To\ grow with time. The velocity perturbation vi is not compared with any unperturbed 
quantity because there is no single unperturbed velocity which could serve as reference at 
all length scales: at scales smaller than the Jeans length, the sound velocity should be used 
whereas one should choose y/AnGpo/k for larger scales. 



A fundamental timescale in the problem is the period of a sound wave, Ts — 2t: /ujs, with 
ujs defined in equation (18). Other important timescales are as follows: 

2.4- 1- Timescales Associated with Background Expansion and with Self-Gravity 
(a) The characteristic time of background expansion is: 



This time is also the time-scale for the adiabatic temperature decrease which follows from 
the expansion (see eq. [7]) and for the inertial change of the perturbed velocity associated 
with the background expansion (first term on the right-hand side of eq. [16]). 

(b) The characteristic growth time of a perturbation by self -gravity, Tg, is 




(22) 



2.4. 



Chciracteristic Timescales 



a 

^ Idl 



(23) 



Tg = (47rGpo) 



1/2 



(24) 



which coincides with cUg^, as defined in equation (18). Prom equation (9), Tg is also the 
timescale for the deceleration of the background expansion. 
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2.4-2. Timescales Associated with Cooling and Conduction 



(c) The characteristic timescale of background cooling because of radiation, Tc, is given by 
(see eq. [7]) 



Tc — I 



-1 



with flc = {l — 1) 



7^^o 



(25) 



For the growth (or otherwise) of the perturbations, one has to consider not only the coohng 
function itself at a given value of Tq and po, but also the differential cooling at neighboring 
values of temperature and density. The precise expression of this fact is obtained through 
the perturbed energy equation (eq. [17]) rewritten in the form: 



(26) 



with the growth rate due to differential cooling, fide, defined as 



(27) 



In the simple case when the perturbation is carried out keeping pressure, density or entropy 
constant, pi/Ti in equation (27) only depends on the background variables. Then, fl'^^ and 
Tc are the characteristic timescales of cooling. 

(d) From equation (26), the timescale of thermal conduction is t^, which was defined in 
equation (19). 



The difi^erent characteristic times have different dependences on k, so the correspond- 
ing processes will become important in different wavelength domains. For instance, the 
conduction growth-rate, t^^^, depends on k as k'^. Hence, conduction dominates at short 
wavelengths. The cooling, background expansion and self-gravity time-scales have no spa- 
tial dependence. Hence, these processes are dominant at long wavelengths. The frequency 
of sound, finally, depends on the first power of k. Hence sound propagation could be the 
dominant effect at intermediate wavelengths. The boundaries between the different regions 
can be easily estimated using the definitions of the characteristic times just given. 



2.5. Solution Method: WKB 



To find the time evolution of the perturbations, the linear ODE system (15)-(17) has 
to be solved. Unfortunately, no simple analytical solution can be found in a general case. 
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The only exception is when all coefficients are time-independent: then, a Fourier analysis 
in time can be carried out and the system reduces to an algebraic eigenvalue problem, from 
which the normal modes of the problem can be obtained. This method was used in many of 
the early papers on thermal instability, like Field (1965). If the background is expanding, 
or is undergoing net cooling, though, the coefficients are time-dependent and the system 
must be solved numerically in general. However, a general statement about the stability 
of the physical system can only be given via analytical solutions, even if they are only 
approximations to the actual solutions. In the following, we study the WKB solutions to 
the problem. 

To carry out the WKB analysis, we first transform the system (15)-(17) into a third 
order linear ODE for pi/po- 



at \po 



Po 



0, 



(28) 



where 



cot 



A{t) = ^(Qep 



(4-37)--fic-^^ 



cp 



(29) 



B{t) = ujI-u' + 2--2- 



(2-37)--f]e-f^cp- 



(30) 



C{t) 



+ ^cp + r~' + (37 + 1) 



(31) 



The WKB analysis requires a clear separation of time scales, in the sense that there is a 
group of processes, which we will call dominant, acting on a time-scale which is much shorter 
than the shortest time-scale of the remaining processes (henceforth called non-dominant pro- 
cesses). In addition, for WKB to be applicable, the time-scales of evolution of the background 
(net cooling, expansion) must be in the non-dominant group. Call To{t) the shortest of the 
characteristic times associated with the dominant processes and Ti{t) the shortest of the 
characteristic times associated with the non-dominant processes, with To{t) -C Ti{t). Call 
coo = TQ^{i) and Ui = rf ^(t), where t is a given time, arbitrarily chosen. The standard 
WKB perturbation method (e.g.. Bender & Orszag 1978; Carrier, Krook & Pearson 1983) 
provides the following formal solution for equation (28): 



— W = — (^o) exp 
Po Po 



E^O->n(t') 



.n=0 



dt' 



(32) 



where the ■0n(^) ^-re unknown functions to be solved for and t^ is the initial time. The growth 
rate of pi/po is thus given by the term in brackets within the integrand. If we substitute 
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equation (32) into (28), an equation can be obtained for each order in cuq- To calculate it, 
we split the coefficients A(t), B{t) and C{t) into components of order uJq, for i = through 
3: 

(A{t) = As{t) + A2{t) + Ai{t) + Ao{t) , 

i B(t) = B2(t) + B^(t) + Bo(t) , (33) 

[c(t) = C,(t) + Co(t). 

The general solution for the functions ipn for arbitrary order n > can be obtained by 
solving the following algebraic equation: 



(Jo 



3— n 



n n—j n— 1 \ 

n—j—m 

vj=0 m=0 j=0 / 

+B2a;(;->, + 5ia;o'->,_i + 5oa;o^->„_2 + Ag-n = . (34) 

In this expression, any function with a negative subindex must be considered equal to zero. 
Equation (34) gives as a function of the previous orders {j < n). The equation for the 
zero-order solution, in particular, is as follows: 

{oooiJof + Ci (cuo^o)^ + B2 ojo^o + A3 = . (35) 

In the special case that propagation of sound, conduction, and differential cooling are the 
dominant processes, equation (35) reduces to the dispersion relation of Field (1965). 

A few properties of the WKB solution are as follows: 

(1) The asymptotic series (32) is convergent if 

k">n+i(i)|«ko-"V'n(i)|. (36) 

These conditions are satisfied if To{t) -C Ti{t) as can be inferred from the relation a;Q~"t/'„(t) ~ 
0[ljo (o^i/tuo)"]) where O indicates the order of magnitude in power of cui/cuq . This relation 
has been obtained using equation (34) and the equations for the background evolution. 
However, the first terms of this series can sometimes give an accurate approximation to the 
exact solution even when To{t) and Ti{t) are of the same order with To{t) < Ti{t). 

(2) The WKB expansion is a global perturbation method, that is, the WKB solution up 
to a certain order is good for all times, not just for times much smaller than the smallest 
background characteristic time (as would be the case for a Taylor expansion). 
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(3) We have assumed a uniform background. In spite of this, the criteria we are going to 
obtain are apphcable locally in a non-uniform background for a Lagrangian fluid element if 

(37) 

where A is the perturbation wavelength and Lq is the macroscopic variation length of the 
unperturbed quantities. That is, the effect of the non-uniform background is less important 
than the effect of the non-dominant processes if condition (37) holds. An additional condition 
to apply the stability criteria in a non-uniform background is: the local uniform expansion 
for the Lagrangian fluid element 

• Vo= - (38) 

should be much larger than the components of the pure deformation tensor and the curl of 
the unperturbed velocity. 

In the following sections, we consider various possible cases of dominant and non- 
dominant processes (see Table 1). In section 3, intermediate scales, for which sound is 
one of the dominant processes, are considered. In section 4, extreme scales, for which sound 
is a non- dominant process, are considered. This includes small scales, for which conduction 
is dominant, as well as long scales, for which self-gravity, cooling or expansion are dominant. 



3. THE INTERMEDIATE WAVELENGTH RANGE: SOUND AS A 

DOMINANT PROCESS 

At intermediate wavelengths, the propagation of sound can be the dominant process, 
either alone or at the same level as one of the other physical processes. In the first subsection 

Table 1: Classification of sections depending on the character of the physical processes 



Section 


Sound 


Self-Gravity Backgr. expansion 


Cooling 


Conduction 


3.1 


D 


ND ND 


ND 


ND 


3.2 


D 


D 


ND 


ND 


4.1 


ND 






D 


4.2 


ND 




D 





Note. — The physical process can be dominant (D) or non-dominant (ND). The hyphen in section 3.2 
indicates that background expansion must be excluded in that section. 
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(sec. 3.1), we consider the first possibility, whereas in the following subsection we study 
the case in which sound and self-gravity are simultaneously dominant in the absence of 
background expansion (sec. 3.2). The readers who are mainly interested in the thermal 
stability of the hot interstellar medium, where self-gravity is negligible, can skip sections 3.2 
and 4 and go directly to section 5 after reading section 3.1. 



3.1. The Sound Domain 

When the sound frequency is much larger than the other characteristic frequencies, there 
are three WKB solutions: two sound-waves and a condensation mode, with growth rates for 
Pi/ Po which we label %l)s± and ■j/'ci respectively. The solutions to zero WKB-order, V-'osi and 
'(/'oc, are obtained by calculating the three roots of the polynomial (35). The coefficients Ai, 
Bi and Cj for this case are given in Appendix A.l, expressions (Al). The result is: 

<^oV'os± = ±i<^s , (39) 

and 

^oc = 0. (40) 

To zero order, therefore, we have the classical sound waves (39) and a non-evolving mode 
(40) which corresponds to any perturbation with = 0, pi = and pi/po 7^ 0. The first 
order corrections, which we label ipis± and ipic-, can be calculated using the general WKB 
solution (34) for n = 1 (and the coefficients Ai, Bi and Cj of eq. [A2]). The results are given 
in the following separately for the condensation mode and for the sound waves. 



3.1.1. The Condensation Mode 



The first order correction for the condensation mode is 



(7 - 1) 

7 7^ 



df [f 



Po 0? 



(41) 



Together with ipQc = 0, the foregoing indicates that this mode does not propagate. Conduc- 
tion has a stabilizing effect, whereas the effect of cooling depends on the sign of [d{C/T)/dT]p. 
The reason for the appearance in equation (41) of the partial derivative oi C/T at constant 
pressure is that pressure remains almost spatially uniform in the condensation mode. Indeed, 
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the solution for pressure and velocity is: 

Vl ^ ^tlcPl ^ nf^\ Pi 

Cs (^s Po V^o J Po ' 

(42) 

Pi _ _ V^ic + 2(a/a)V'ic + -01c - uj^g pi nf^\ £l 
Po - ^ Po ^ ^ [co^J Po ' 

where O indicates the order of magnitude in power of cui/cuo, with cui in this case given by 
the shortest of the conduction and cooling times appearing in equation (41). Hence, the 
relative pressure perturbation obtained from the linear system is zero down to second WKB 
order. 

The second order WKB correction for the growth rate of pi/po for the condensation 
mode is equal to zero as can be shown using equation (34) for n = 2: the following non-zero 
correction appears at third WKB order. 



3.1.2. The Sound Modes 



The first order WKB correction for the sound modes is: 



(7 



P I 7 



n 



7 



df 



Xo k"^ 



+ — — 
Po a 











1) 


1} 



(87 — 5) a 
4 a 



(43) 



where 5* stands for the specific entropy. Since -^os has no real part, the change in amplitude 
is given by ipig. From equation (43) we see that a net background radiative cooling (£0 > 0) 
has a destabilizing effect. Expansion tends to destabilize these modes if 7 > 5/3. Conduction 
has an unconditionally stabilizing effect for the sound mode. In contrast to the condensation 
mode (see eq. [41]), now the T-derivative of C that appears in the growth rate is at constant 
entropy. This is easily understandable, since entropy is changed by conduction and cooling 
only, which are non-dominant processes in the present case. So, these modes evolve essentially 
at constant entropy. Self-gravity does not appear in this first order correction (assuming 
that all the non-dominant processes have a characteristic time of the same order). Putting 
together the zeroth- and first order expression, we can write the WKB solution for the density 
perturbation in the sound mode as: 
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The solution for the pressure and density up to first WKB order for the sound waves is: 
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To zero-order, this reduces to the classical pi/po = {^h){Pi/Po) = T^i /cs of the elementary 
sound waves. 



The solutions of the previous subsections generalize those obtained by previous authors 
for particular cases. Thus, in the absence of background expansion and with no net back- 
ground cooling (i.e., when there is thermal equilibrium in the unperturbed state), our WKB 
solutions (39), (41) and (43) reduce to the solutions of Field (1965). In the absence of 
thermal conduction, for 7 = 5/3 and assuming that C/p only depends on T, our solution for 
the condensation mode reduces to that of Mathews & Brcgman (1978). In the absence of 
thermal conduction and for 7 = 5/3, our solutions reduce to those of Balbus (1986). 

For the temperatures and densities for which cooling is destabilizing. Field (1965) 
found that for each mode, sound and condensation, there exists a critical length, which we 
will call henceforth classical Field length, defined as the maximum wavelength for which 
conduction can suppress the instability. Field considered a static background in thermal 
equilibrium. For a background undergoing cooling and expansion, wc can then now define a 
generalized Field length (and the corresponding generalized Field wavenumber, kcrit) for each 
mode as the wavelength (or wavenumber) separating stability domains. The generalized 
Field wavenumber for the condensation mode is thus given by: 



This expression has been obtained setting equation (41) equal to zero. Note that cooling 
must have a destabilizing effect {[d{C/T) / dT\p < 0) for the existence of this critical length. 
There is also a generalized Field length for the sound waves, given as the positive real solution 
of equation -^ig = which usually is larger than the one for the condensation mode. 



3.1.3. Link with the Previous Literature: the Generalized Field Length 




(47) 



3.1.4. Comparison of WKB and Exact Solutions for a Particular Case 



In this section, we consider an idealized case to compare the WKB solutions and the nu- 
merical solutions of the system (15)-(17) (and leave for section 5 the apphcation to a realistic 
ISM or IGM). We consider a completely ionized plasma of pure hydrogen with temperature 
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above lO^K. For those temperatures, cooling is almost only due to bremsstrahlung because 
the line and recombination emission due to elements heavier than hydrogen is not present. 
Conduction, in turn, is due almost exclusively to electrons. Thus, L and x can be taken in 
the form (Rybicki & Lightman 1979; Spitzer 1962): 

C = lO-^njyT^/2 erg gr-^s"^ , (48) 
1 84 X lfl~^ 

X = — ^ 1 T^''^ erg s-^K-^cm-i , with log Ac = 29.7 + log (T x 10-Vne^/^)(,49) 

log Ac ^ ' 

where the averaged Gaunt factor has been chosen equal to 1.2, nn is the proton number 
density in cm~^, is the electron number density in cm~^, and T is measured in K. In this 
subsection, the Coulomb logarithm log Ac is set equal to 29.7. The cooling function (48) 
has a destabilizing effect on the condensation mode (see eq. [41]) and on the sound waves 
(see eq. [43]). We assume a non-expanding background, so that po remains constant in time. 
The evolution of the background temperature can be easily obtained from equations (7) and 
(48). 

Using as time unit the period of a sound wave at time t = 0, Tsq, and equations (48) 
and (49) for cooling and conduction, it turns out that the coefficients of the relative pertur- 
bations in the system (15)-(17) depend on k and po through the single parameter p^jk only 
(equivalently, A uh). This allows a large simplification in the presentation of the results in 
the following. 

(a) Critical wavelengths. The critical wavelength for the condensation mode (i.e., the 
generahzed Field length) is obtained from equation (47) , whereas for the sound waves it can 
be calculated by setting the growth rate (43) equal to zero. Both are represented in Figure 
1 as implicit functions of Anjy and Tq. The critical wavelength for the condensation mode 
(solid line), is smaller than the critical wavelength for the sound waves (dashed line). Figure 
1 also shows the region (drawn as a shaded band) where r^~^ < ujg and uj^ > t~c^ ■ The 
sound domain is contained within that band, but, as seen just below, the WKB solutions of 
the previous subsections are good approximations also close to the boundaries of the region. 
Well above the shaded band, cooling acts on the fastest timescale; well below it, conduction 
is dominant. 

(b) Comparison of WKB- and exact solutions. We consider both the condensation 
mode and the sound waves for a case defined by Xhh = 10^^ cm~^ and Tq = 7.3 x 10^ K at 
t — Q (Fig. 2). With this choice, both modes have initially a wavelength below the critical 
value. Figures 2a and b show the time evolution of the background temperature and of 
|^cp|/<^s and jiOg. The other four panels in the figure compare the WKB solution with the 
numerical solution of the linear system (15)-(17) for an initial amplitude of the perturbation 
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Fig. 1. — Critical wavelengths and boundaries of the sound domain. The solid line indicates 
the critical wavelength for the condensation mode. The dashed line corresponds to the 
critical wavelength for the sound waves. The two dash-dotted lines are the points where 
— cUg (upper straight) and ^ = Us (lower straight): the sound domain, therefore, is 
contained within the shaded band. 
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Pi/po — 0.01. Panels (c), (d) and (e) show the time evolution of the condensation mode 
(solid: WKB approximation; dashed: numerical solution) and the relative error between 
them (dotted line). The WKB solution turns out to be an excellent approximation to the 
exact solution: both curves are almost indistinguishable except when \Vtcp\ ~ cu^. In fact, 
the relative error of the WKB solution for pi/po is still below 0.01 when \VLcp\/^s reaches 
the value 0.1. Initially, the condensation mode has a wavelength below the critical value, 
but the latter decreases as the system cools (see Fig. 1) and their ratio grows above unity 
for t/TsQ — 12. Therefore, pi/po first decreases, then reaches a minimum at that time, and 
grows thereafter. Figure 2f shows the WKB and numerical solutions for the sound wave. 
Again, the agreement between both solutions is excellent. The wavelength ratio follows a 
similar evolution as for the condensation mode, and so does the perturbation amplitude. 



3.2. Solutions on Scales around the Jeans Length: the Acoustic— Gravity 

Domain without Background Expansion 

Going to spatial scales larger than in the foregoing section, there is interesting new 
physics when the wavelength approaches the Jeans length, Lj — I'KCsjojg , that is, as soon 
as LOg is no longer small compared to 0;^. In this subsection we study the regime in which 
either sound or gravity (or both simultaneously) are dominant and the background is slowly 
cooling or heating but with zero expansion rate (if d 7^ 0, the expansion would become as 
dominant as self-gravity itself after a time of order Tg at the latest, see eq. [9], and the WKB 
model could not be used). 

Consider, then, the case in which \^dc\^ and are all much smaller than both lo^ 
and LOg. There are then three WKB solutions: two acoustic-gravity modes and a condensa- 
tion mode. We use the symbols ^ag-t and for their respective growth rates for pi/po- As in 
the previous case, the zero-order solutions are obtained by calculating the three roots of the 
polynomial (35), but now the coefficients Aj, Bi and Cj have the values given in Appendix 
A. 2, expressions (A5). The result is: 

^oV'Oag± = ±i^u1-L0l (50) 

and 

V'oc = 0. (51) 

To zero order, therefore, we have a couple of acoustic-gravity modes (eq. [50]) like those 
obtained by Jeans (1902) and a non-evolving mode (eq. [51]). The latter corresponds to 
any perturbation with f i = and Pi/po = '-f{ujg/ujf)pi/pQ, i.e., equilibrium between pressure 
gradient and self-gravity. If ujg < ujs the acoustic-gravity modes are waves. They have 
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Fig. 2. — Comparison of the WKB and the numerical solutions for the condensation and 
sound modes, (a) Evolution of the background temperature, (b) Solid line: time evolution 
of |r^cp|/'^s- Dashed line: evolution of /iOg. (c), (d), and (e) show the evolution of the 
condensation mode: the solid hne is the WKB solution and the dashed line is the numerical 
solution (in panel d, they coincide within the resolution of the figure). Panel (f) shows the 
evolution of the density perturbation for a sound wave (again, WKB and numerical solution 
coincide in the figure). Panels (c), (d) and (f) also show, as a dotted line, the relative error 
between the numerical and the WKB solutions. 
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a frequency smaller than a pure sound wave with the same wavelength because the self- 
gravity force acts against the pressure gradient force. If uig > Ug these modes grow or decay 
exponentially. In the growing mode, in particular, the perturbation collapses gravitationally 
and the pressure gradient is not large enough to prevent it. The introduction of cooling and 
thermal conduction as non-dominant processes makes the amplitude of these modes change. 
These modifications are given by the WKB corrections of order one and greater. 

To calculate the first order corrections for the growth rates of Pi/po for those modes, 
namely ipiag^ and ipic^ we turn again to the general WKB solution (cq. [34]), using n = 1 
and the coefficients given in expressions (A6) of Appendix A. 2. We show and discuss them 
in the following: 



3.2.1. Condensation Mode 
The first order WKB correction for the condensation mode is given by 



d_ fC 
df \T 



4(1) +i£^!(i-4 




(52) 

This mode reduces to the sound domain condensation mode (41) when ujg <^ Ug. We 
introduce the quantity 



F = ^ 74- . (53) 

which is proportional to the perturbation of the net volume force (pressure gradient minus 
self-gravity). Using this quantity, the relative perturbations can be written 

v_.^i^Pi^^(^\pi^ F = -^tl^Pl^o(i)^, (54) 

Cs UJs Po \i^sj Po Po \OJsJ Po 

where uji is the inverse time-scale of the non-dominant processes (conduction and cooling). 
This is similar to equation (42) but now the sound-gravity modes, that propagate or grow 
much faster than the condensation mode, lead to the cancellation of the gravity force by the 
pressure gradient instead of to a condition of uniform pressure throughout. 

The first order correction (52) shows three interesting features. First, it diverges when 
ujg — cug, i.e., there is a turning-point in the WKB sense. The turning point can easily be 
reached in a cooling and non-expanding medium for wavelengths for which, initially, cus > t^g'- 
the former frequency decreases with time whereas the latter stays constant. As known from 
the WKB theory, when a solution goes through a turning point, a mixing of modes takes 
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place. Second, conduction has a destabilizing effect in the range 70;^ > u>g > Ug . This can 
be explained using the equation of evolution for the perturbed entropy, 5*1, and the relation 
between Si and Ti/Tq following the condition of permanent equilibrium between gravity and 
pressure gradient. Third, the combination of partial derivatives of C is quite different to the 
corresponding expression in equation (41). In the following we explain the origin of the last 
feature. 

To explain the origin of the cooling terms in equation (52) we can use the results of 
Balbus (1986). Using Lagrangian perturbations in the entropy equation, he derived an 
instability criterion for the condensation mode in the absence of conduction in the particular 
case in which the perturbation of an arbitrary thermodynamic variable A is kept equal to 
zero during the evolution of the condensation mode. The resulting criterion is obtained from 
the following equation: 

1 dAS 



AS dt 



\- 




_dS 


m 



(55) 



A 

where A indicates Lagrangian perturbation. Balbus (1986) applies this equation to the 
cases A = p and A = p. In the present section, in turn, the physical quantity whose 
perturbation remains equal to zero is F (eq. [54]); on the other hand, Lagrangian and 
Eulerian perturbations are indistinguishable since the background is uniform (e.g., Shapiro 
& Teukolsky 1983). It can then be shown that the cooling terms in equation (52) can be 
derived from the right-hand side of equation (55), using the relationship between the growth 
rates of Si and pi/po- 

Critical lengths are computed solving the equation ipic = 0. This equation can be 
expressed as: 




X ^ T 1 O \ \ 2 1 r 



Jlp 7 j Lij 7 Qp 

where \crit is the critical wavelength, r ~^ is evaluated for A = Lj, and 



(56) 



Q. = (T-l)^fi) . (57) 



n \dT J 



p 



i.e., minus the growth rate due to isochoric cooling in the classical theory of Field (1965). 
At the Jeans length, it is natural to expect cooling to be much faster than conduction 
<^ Of). Assuming this, the solutions of equation (56) are given by: 

^critX = 77 -^J and AcHt^2 = ~7T^-^J • (^S) 

For the existence of a critical length, the right-hand side of the corresponding expression must 
be positive. The first critical length coincides with that of equation (47), which corresponds 
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to a condensation mode in the sound dominated case. The second critical length is new, and 
has a value that is near Lj. If both critical lengths exist and \criti,2/Lj < 1, the condensation 
mode will be unstable for wavelengths A such that Xcriti < A < Xcrm- 



3.2.2. Acoustic-Gravity Modes 



The first order WKB correction for the acoustic-gravity modes is 
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(59) 



This term coincides with the first-order correction for the sound waves given in equation 
(43) for zero background expansion but divided by 1 — ujg/ujg. For Ug < Ug, this is not 
surprising because the acoustic-gravity waves behave in that case like sound waves of the 
same wavelength but with a smaller frequency. Similarly as for the sound domain, for this 
mode there is only one critical wavelength, obtained solving ipiag — 0. 



4. THE SHORT- AND LONG-WAVELENGTH RANGES 

For completeness (and later use in section 5), we briefiy discuss in this section two 
results for perturbations in the conduction and cooling domains. 



4.1. Small Spatial Scale: Heat Conduction Domain 

At small enough spatial scale (a) thermal conduction is the dominant process and (b) 
uj~^ is much smaller than Tg, Tg and Tc, so that ^ cj^T^ ^ min(rg, t^., r^). Two simple 
consequences of the foregoing are: (1) conduction can eliminate any temperature gradients 
well before the sound waves can cross a wavelength and (2) the evolution of the dominant 
processes (conduction and sound propagation) occurs without important evolution of the 
background through coohng and expansion. The solutions to zero and first WKB order 
must thus coincide with the solutions obtained by Field (1965), viz. a condensation mode at 
zero order and quasi-isothermal sound waves at first order. The condensation mode evolves 
on the timescale of the thermal conduction, r^. Using this fact in equations (15)-(17) one 
obtains vi/cg ~ 0{uJs/t~'^) Pi/po , pi/po ~ 0{uj^/t~'^) pi/po and 

d^P_.\_]^P_^ 
dt \poJ T^Po 
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The condensation mode at small length-scales is therefore completely stabihzed by conduc- 
tion. As a consequence of the smallncss of pi/ po, to obtain the WKB solutions for this mode 
it is necessary to use a third order differential equation for pi/po instead of the differential 
equation (28) for pi/ pQ. 



At large enough wavelength, cooling, expansion and/or self-gravity dominate, so that 
the following ordering of timescales applies: min (r^, Te, Tc) ^ cj"^ ^ r^. Whenever cooling 
or expansion dominate, no WKB solutions can be calculated, since in that case background 
and perturbation evolve on the same or comparable timescales. However, a zero order 
solution can be obtained by direct inspection of equations (15) through (17) when there is a 
single dominant process. For instance, if cooling is the single dominant process, we expect 
the zero order solution to have a growth rate of order = max (r~^, VLdc) for ^ relative 
perturbations pi/ po, Vi/cs and Pi/po- Inserting that timescale into equations (15)-(17) one 
finds: Vi/cg ~ 0{u!s/t~1) pi/po and pi/po ~ ^i^s/'^mc) Pi/Po- This perturbation evolves at 
almost constant density because the propagation of sound is much slower than the growth 
by cooling. The growth rate at zero order can be obtained from the previous considerations 
and equations (17) and (20): 



Balbus (1986) obtained this growth rate for the perturbed entropy using equation (55) with 



5. GENERALIZED FIELD LENGTH IN AN ASTROPHYSICAL COOLING 

MEDIUM 



We consider the thermal stability of an optically thin astrophysical medium at a temper- 
ature in the range 10^ — lO^K which is undergoing net cooling. The medium is assumed to be 
uniform (or only weakly non- uniform in the sense of eq. [37]). Old SNRs and super- bubbles 
are examples of such a medium. The generalized Field length, which can be computed using 
the results of section 3.1.3, separates the stable regions from the unstable ones in wavenum- 
ber Fourier space. Self-gravity can be ignored because for the considered medium the Jeans 
length is usually much larger than the Field one. Assuming no heating and permanent ion- 
ization equilibrium, the cooling function only depends on the instantaneous properties of the 



4.2. Large Spatial Scale: Cooling Domain 




(61) 



A^p. 
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medium and can be expressed as follows: 

>C=(^)pA(r), (62) 

where rie is the electron number density, uh is the total number density of hydrogen (atoms 
and ions) and A is a function that only depends on temperature. We use the function A(T) 
provided by J. C. Raymond (2001, private communication) for standard abundances and 
uh — 1 cm~^ (see Fig. 3), which is an improved version of the cooling function of Raymond, 
Cox & Smith (1976). For T < 3 x lO^'K, the function A(r) has a very weak dependence 
on Uh, which can be safely neglected in the present context. In the considered temperature 
range, 7 = 5/3 and the thermal conduction coefficient, x, is given by equation (49). A 
number of authors have suggested a substantial reduction of x in the Intcrgalactic Medium 
due to, for instance, kinetic instabilities driven by temperature gradients (e.g., Pistinner, 
Levinson & Eichler 1996), a tangled magnetic field, which makes the effective mean free 
path of electrons smaller (e.g., Rosner & Tucker 1989) and may cause magnetic mirroring 
(e.g., Chandran et al. 1999). Thus, as is usual, we include a constant factor ^, with value 
between and 1, multiplying the conduction coefficient 



5.1. Unmagnetized Background 

The generalized Field length, Xcrit, can be calculated inserting the cooling function (62) 
into the growth rates (41) and (43) for the condensation and sound modes, respectively, and 
setting them equal to zero. The result is: 

-1 

(63) 

with a = 2 for the condensation mode and a = 3/4 for the sound waves. In equation 
(63), X = ne/fiH and all quantities are evaluated in the background. If the logarithmic 
temperature gradient of A is positive and steep enough, the right-hand side of equation (63) 
becomes negative and there is no critical length: the corresponding mode is stable at that 
temperature for all wavelengths. 

It is instructive to compare equation (63) with the classical Field length (sec. 3.1.3), 
i.e., the critical length if the background were in thermal equilibrium through the action 
of a heating which exactly compensates the cooling. Since we are interested primarily in 
the destabilizing effect of the cooling, we assume a density- and temperature-independent 
heating. The classical Field length is also given by equation (63) but with a = 1 for the 
condensation mode and a — —3/2 for the sound waves. From these values of a, comparing 
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Fig. 3. — Cooling function A(r) of J. C. Raymond (2001, private communication) for stan- 
dard abundances. 
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them with the values obtained above (see Fig. 4), we gather that there are more (and wider) 
unstable temperature ranges in the case with net background cooling than in the classical 
Field problem. In fact, the difference is quite striking for the sound waves, which change 
from stable to unstable almost everywhere if one suppresses the heating providing thermal 
equilibrium. 

Figure 5 shows the instability regions and the critical lengths in the two-dimensional 
parameter space (T, A njj x^^"^ ■C^^^^)- The upper diagram corresponds to the sound waves, the 
lower one to the condensation mode. In the dark-grey shaded region, both the classical Field 
problem and the case with net background cooling arc unstable, whereas in the light-grey 
shaded region only the latter is unstable. As expected, the instability region for the sound 
waves is much larger for the case with net background cooling. For the condensation mode, 
both regions are comparable, but, again, the range of unstable temperatures is smaller for the 
classical Field problem. In fact, the ranges of instability for sound waves and condensation 
mode arc similar when there is net background cooling, which is in strong contrast to the 
situation in the classical Field problem. The upper left corner of each diagram corresponds 
to a region where cooling is dominant (i.e., it occurs on the fastest time-scale of the problem). 
Hence, there are no properly defined sound waves in that region (hatched area in the upper 
diagram) . 

The critical length expression (63) is valid only when sound is the dominant process, 
which happens within the region delimited by the curves ujg = r^^ (dotted line) and = 
(dash-dotted line) in the diagram. However, we know that the region below the = - 
line is stable, since conduction is dominant. Similarly, in the cooling-dominated region (i.e., 
well above the dash-dotted line), where conduction is negligible, we know the instability 
properties of the condensation mode both for the classical Field problem and for the case 
with net background cooling using an isochoric criterion as explained in section 4.2. The 
remaining unknown region is, thus, a (possibly narrow) band around the cUg = - line. 
We have tentatively indicated this band in the figure as a cross-hatched region. As a final 
remark, note that the use of the parameter A tih x 

1/2 ^-1/2 ^J^g g 

gure is meaningful because 

we have ignored the weak dependence of the Coulomb logarithm on (see eq. [49]). For 
the figure, rie has been set equal to lcm~^ in A^. 



5.2. Inclusion of a Background Magnetic Field without Dynamical Effects 

The astrophysical media considered in the previous subsection (like the hot interstellar 
medium or the intergalactic medium) are threaded by magnetic fields. Usually, these fields 
are dynamically unimportant, that is, in these media j3 » 1, where j3 is the ratio of 
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10" 10" 10" 10^ 10 

T (K) 



Fig. 4. — The solid line shows (c? log A)/ (o? log T) for the function A(T) represented in Figure 
3. The horizontal hnes give the upper hmits of existence of a critical length for the conden- 
sation mode (dashed) and for the sound waves (dotted) , whereby in both cases the lower line 
corresponds to the Field problem and the upper one to the case with net background cooling. 
Wherever the sohd line falls onto the shaded bands, the system with net background cooling 
is unstable whereas the results for thermal equilibrium predict stability. 
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Fig. 5. — Critical wavelengths and regions of instability for the sound waves (upper panel) 
and the condensation mode (lower panel) for the application discussed in section 5. The 
critical wavelength \crit is plotted as a solid line for the case with net background cooling 
and as a dashed line for the classical Field problem. The dark-grey shaded area is the region 
where both problems yield instability; in the light-grey shaded region only the case with net 
background cooling is unstable. In the upper left corner of the upper panel (hatched area in 
the figure) the sound waves are not properly defined so that no statement concerning their 
stabihty is possible there. Two further lines depicted in the figure are the locus of points 
where a;^ = t~^, with ^ = 1 and x = 1 (dotted) and those where oUg — t~^, with ^ = 1 and 
X = 1 (dash-dotted line). The stability of the points in the figure in a narrow band around 
the dash-dotted line (indicated with cross hatching) is uncertain, as explained in the text. 



-29- 



thermal and magnetic pressures. However, the effect of this weak magnetic field on thermal 
conduction is very strong. Therefore, we include such a weak magnetic field to complete the 
discussion of the previous paragraph. In a medium threaded by a magnetic field, the growth 
rate of pi/po foi' the condensation mode in the limit /? >> 1 is given by (A. J. Gomez- Pelaez 
& F. Moreno-Insertis, in preparation): 
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where 6 is the angle between the background magnetic field, Bq, and the wavenumber vector 
k, the conductivity parallel to the magnetic field, X|h given by equation (49), and the 
perpendicular conductivity, is given by (Spitzer 1962): 
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in which logA^; is given by equation (49), rij is the number density of ions in cm~^, B^q 
is the magnetic field in //G, and Tq is the temperature measured in lO^K . Expression (65) 
is only applicable when the temperature of ions and electrons is the same and the Larmor 
radius of the ions is much smaller than their mean free path, which is always the case in the 
media considered in this section. The corresponding expression for the sound waves is: 
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In the limit P » 1, the only effect of the magnetic field is to make the thermal 
conduction depend on the spatial direction. The growth rate of both modes, and therefore 
their associated generahzed Field lengths (Xcrit): depend on 6. Setting both, t/jic (eq. [64]) 
and ipis± (eq. [66]), equal to zero, one obtains 



cos' e + — sin' 9 
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(67) 



where Xcrit\\, is given by equation (63). We will also use the symbol Xcrit± for Xcrit{9 — 7r/2). 

Figure 6 shows the generalized (sohd fine) and classical (dashed fine) Field length and 
the instability regions for the condensation and sound modes with 9 — 7r/2. The upper 
diagram corresponds to the sound waves, the lower one to the condensation mode. Here 
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again, dark-grey shading indicates that both the thermal equihbrium problem and the case 
with net background cooling are unstable; in the light-grey shaded area only the latter is 
unstable. Note that now the combination XB^qx'^^^ appears in ordinates, since Xcriti. does 
not depend on Uh but on -B^g*- From the figure (compare with Figure 5), it is apparent 
that Xcriti. << ^crit\\ and that Xcriti. has a weaker dependence on temperature than Xcrit\\- 
Both features are a direct consequence of the conductivity ratio (65). Hence, if there is a 
background magnetic field, the condensations, in the linear stage, will be filaments directed 
along the magnetic field lines with a length of the same order as Xcrit\\ and a diameter of the 
same order as Xcriti_- 

Typical values of the critical length in the hot interstellar medium (as inside SNRs 
and super bubbles) follow from the physical parameters in them like: uh — 10"^ cm~^, 
T = 10^ K, and B^a = 1 (e.g., McKee 1995). The associated critical lengths are: (1) 
^crit\\ — 3.4 kpc and Xcrit± = 9 x 10~^ pc for the condensation mode and (2) Xcrit\\ — 6 kpc 
and Xcriti. = 1-5 x 10~^ pc for the sound waves. Note that Xcrit\\ is larger or of the same 
order than the size of the considered systems whereas Xcrit± is well below it. 

6. SUMMARY AND DISCUSSION 

We have carried out a systematic analysis of the thermal stability of a uniform medium 
which, in the unperturbed state, is undergoing net cooling and expansion, so that the prob- 
lem is not amenable to a full Fourier normal mode analysis. Thermal conduction and self- 
gravity have also been included as fundamental ingredients. The small-perturbation problem 
yields a system of ordinary differential equations with three independent solutions. In many 
cases, two of them have the character of oscillatory modes (typically sound waves, or, more 
generally, acoustic-gravity modes) and the third one is a non-oscillatory (or condensation) 
solution. The wavelength of the perturbation determines the positions of the different physi- 
cal processes in the hierarchy of timescales. At small enough wavelength, thermal conduction 
damps the condensation mode and makes sound waves isothermal. At large wavelength, one 
or a few of the following processes are dominant: cooling, self-gravity and expansion. For 
example, if cooling is the only dominant process, the perturbation evolving on the fastest 
timescale is isochoric with no soimd waves present that could even out the pressure gradi- 
ents. At intermediate wavelength, the propagation of sound is the fastest process, so the 
evolution of the condensation mode is isobaric. In the paper, solutions have been obtained 
in a number of instances using the WKB method. WKB solutions can be found only if two 
conditions are fulfilled: (1) there is a clear separation of timescales between fast and slow 
physical processes in the system, and (2) the background expansion and cooling belong to 
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Fig. 6. — Perpendicular Field length, Xcrit±-, a-nd regions of instability for a magnetized 
medium with /3 ^ 1 for the sound waves (upper panel) and the condensation mode (lower 
panel). Xcrit± is plotted as a solid line for the case with net background cooling and as a 
dashed line for the classical Field problem. In the dark-grey shaded areas both the thermal 
equilibrium problem and the case with net background cooling are unstable; in the light-grey 
shaded area only the latter is unstable. 
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the slow category. 

Of special interest in this paper are those solutions whose growth rate is of the same 
order as the rate of change of the unperturbed background. In this case, the change in time 
of the critical wavelengths (the generalized Field length, Xcrit, and the Jeans length, Lj) may 
make the perturbation go through totally different stability regimes during their evolution. 
For example, the solutions obtained in section 3.2 for a medium undergoing net cooling 
include a condensation mode initially in the stable range whose wavelength becomes larger 
than the Field length after some time, thus becoming thermally unstable. Later on, the Jeans 
length decreases sufficiently so that the perturbation also becomes gravitationally unstable. 
A sound wave may undergo a similar process of change of stability, first being in a stable 
regime with slowly decreasing amphtude, then becoming thermally unstable and finally, 
when the self-gravity frequency becomes comparable with the sound frequency, turning into 
an acoustic-gravity mode which ends up as a gravitationally-dominated perturbation. 

The analysis of the present paper shows (section 5) that in a medium undergoing a realis- 
tic net cooling the sound waves are thermally unstable for a comparable range of wavelengths 
and temperatures as the condensation mode. This may come as a surprise: in contrast, and 
by way of example, if the thermal equilibrium is reached by balancing the cooling with a 
constant heating, the sound waves are stable for almost all temperatures in the range 10^ 
- 10* whereas the condensation mode continues being unstable. On the other hand, the 
inclusion of a weak magnetic field has a strong infiuence on the thermal conduction. For 
astrophysical environments such as the hot ISM and the hot IGM, the generalized Field 
length along the magnetic field is many orders of magnitude larger than in the perpendic- 
ular direction. Therefore, the unstable condensations will be filaments directed along the 
magnetic field lines. 

For the description of the astrophysical systems in which thermal instabilities appear, 
it is important to study the non-linear phase of both condensation mode and sound waves. 
The unstable sound waves develop pairs of shock fronts and rarefaction waves, as known 
from elementary hydrodynamics. The condensation mode, instead, just continues growing 
without propagating. The detailed physics of the nonlinear phase has to be calculated, in 
general, using numerical means. Yet, a number of results can be advanced at this stage. For 
instance, in the linear analysis, a perturbation with negative pi/po grows (in absolute value) 
as fast as a perturbation with positive pi/po- This symmetry is broken in the non-linear 
evolution. In a medium undergoing net cooling, for example, a condensation with positive 
Pi/po (negative Ti/Tq) grows faster than one with pi/po < 0. This happens because of 
the sum of the following two processes. First, thermal conduction grows with temperature. 
Therefore, a condensation with positive Ti/Tq is more damped by thermal conduction than 
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a negative one. Second, in the range 10^ — 10^ cooling increases toward lower temperatures 
(Fig. 3). Therefore, a condensation with negative Ti/Tq is more destabilized by cooling than 
a positive one. A related aspect is that the characteristic spatial size (width at half height) of 
a condensation decreases by orders of magnitude during its evolution. This happens because 
the coldest zones of the perturbation cool much faster than the warmer ones. The last aspect 
can be seen in the paper of David, Bregman & Seab (1988), where the non-linear evolution 
of a perturbation of this type is computed numerically. 

In the weakly-magnetized case, at the beginning of the non-linear evolution {j3 ^ 1), the 
compression in the filament proceeds mainly in the directions perpendicular to the magnetic 
field lines. This is because of the reduction of thermal conduction transversely to the field 
lines: a much higher temperature (and thus pressure) gradient can be maintained across than 
along the field fines. The plasma /3 inside the condensation decreases with time because of 
the temperature decrease and the growth of magnetic field due to the compression. In the 
advanced non-linear phase, the perpendicular compression comes to a halt once /3 <^1, since 
then the net force across the field lines (the gradient of the sum of magnetic and thermal 
pressures) is almost decoupled from temperature. After this, the compression is only possible 
along the field lines. As a related example, see David & Bregman (1989), which numerically 
compute the non-linear evolution of a planar condensation taking into account magnetic 
pressure. 

The WKB analysis that we have carried out in this paper assumes a uniform background. 
However, our results are also applicable for a weakly non-uniform background. More in detail, 
our results for the condensation mode in the sound domain are valid if the Brunt- Vaisala 
frequency associated with the background entropy gradient is much smaller than the growth 
rate of the condensation mode. This holds, for example, in the supersonic region of a wind 
and inside old SNRs and superbubbles. In contrast, for the study of the thermal instabifity 
in cooling fiows in clusters of galaxies, the buoyancy of the perturbations must be taken into 
account. Many authors have studied this problem: Mathews & Bregman (1978), Malagoli, 
Rosner & Bodo (1987), Balbus & Soker (1989), Loewenstein (1990) and Balbus (1991). 
Loewenstein (1990) and Balbus (1991), in particular, conclude that a magnetic field, even 
a weak one, can inhibit buoyancy and trigger thermal instability, which otherwise would be 
small due to buoyancy. The extension of the results of the present paper to environments of 
this type requires additional work. 

Financiation through DGES project PB95-0028C of the Spanish Ministry of Education 
is acknowledged. The authors are grateful to Dr J. C. Raymond for providing them with the 
cooling functions used in this paper. 



COEFFICIENTS FOR THE WKB SOLUTION 



A.l. Coefficients for the Sound Domain 

For the sound domain case (see sec. 3.1), the different order components of the coeffi- 
cients A{t), B{t), and C{t), have the value: 



^3 = 0, B2 = ujI, Ci = , 



(Al) 



A2^'^{n,p + T^-') , B, = 0, Co = f^c + ^cp + V' + (37 + l)^, (A2) 



Bo = 2-- u:l - 2- 
a ^ a 



(2-37)--l],-l]ep-r^-' 
a 



(A3) 
(A4) 



A. 2. Coefficients for the Sound and Self-Gravity Domain without Expansion 

When sound and self-gravity dominate and there is no expansion (see sec. 3.2), the 
different order components of the coefficients A{t), B{t), and C{t), have the value: 



^3 = 0, B2 



2 2 



Ci = 0, 



^2 = Y (^cp + T^"') - uj] + l^ep + r^"'] , 5i = , Co = ^c + ^cp + r^"' • (A6) 



(A5) 
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